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ABSTRACT 


A theory is presented for computing the far field beam 
patterns from distributed random noise sources. The theore- 
tical model utilizes the Green's Function for the wave equa- 
tion and the space-time autocorrelation function for deter- 
mining the radiation from a randomly vibrating area. The 
actual far field beam pattern of a horn speaker in an anechoic 
chamber was obtained, and also near field measurements were 
taken to obtain the correlation distance and the mean square 
of the particle velocity using the autocorrelation function. 
Finally a computer program was written to evaluate the inte- 
gral wave equation by numerical methods. It was found that 
the critical parameters in the mathematical model were the 
correlation distance and the frequency limits of nite ometeclone 
Small variations in the correlation distance modified greatly 
the width of the predicted beam pattern, while changes in the 
limits of integration had a moderate effect. The Frequency 
Spectrum was oireelinedl in the anechoic chamber and it was used 
to determine the limits of integration of the integral solu- 
tion for the intensity field. 

The predicted beam patterns were in satisfactory agreement 
With the measured far field beam patterns. The differences 
between them are probably due to inaccurates in determining 


the limits of integration and the correlation distance. 
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I. INTRODUCTION 


It has long been recognized that the best means of locat- 
ing a deeply submerged submarine is by detecting the sound 
that it radiates into the water; also all surface ships are 
Subject to attack by passive homing torpedoes and acoustic 
mines, both of which operate using the radiated noise of the 
ship. These facts make the far field radiated signatures a 
very important factor to determine detectability characteris- 
tics and capacity of survival in a hostile environment. 

Free field measurement in the far field involves serious 
difficulties like ambient noise, multipath, convergence zone, 
attenuation, among others, which makes impracticable such 
kind of measurement. Instead a prediction of the Far Field 
uSing Near Field data is suggested. Certain general princi- 
ples are repeatedly invoked in prediction. The first is that 
the vibratory powers of several sources is additive. This is 
based upon the assumption that the signals from any two 
sources are statistically uncorrelated, an assumption which 
is almost invariably found correct. Signal contributions from 
non-overlapping pen tency bands are inherently uncorrelated; 
hence the total power over all frequencies is found by summing 
contributions from various frequency bands. 

several works have been done about prediction of Far 
Field beam pattern from Near Field measurements for discrete 


frequencies (Refs. 1,2,3,4). In general these methods deal 





with the amplitude and phase information taken at the Near 
Field and the final result shows that the Fourier Transform 
of those measurements lead to the beam pattern. Consider 

the random noise as the random summation of sinusoidal waves, 
thus a complete analysis of the frequency bands can be 
achieved using a very narrow filter which center frequency 
can be shifted to cover all the spectrum, or a bank of narrow 
filters can be used. Adding the contribution of each fre- 
quency it should be possible to obtain the final beam pattern, 
but either procedure involves a tedious work, specially for 
broad band signals. There also remains the question of cor- 
relation between different bands. 

The more general problem of predicting the Far Field noise 
beam pattern from a random noise source is of great practical 
importance. The theory of the noise field for systems that 
Satisfy the wave equation has been developed by some inves- 
tigators working in the fields of acoustics (Ref. 5), and 
optics (Ref. 6). Different methods had been suggested in or- 
der to generate a prediction in the far field, one of them 
consists of using correlation functions and the Green's func- 
tion for the wave equation, for general time variations. This 
technique is used by Morse and Ingard (Ref. 7) in their book 
"Theoretical Acoustic," for describing non-periodic vibra- 
tions or random processes. 

This paper discusses one possible method of Near Field 
analysis using the spatial correlation function and the 


Green's function which enables differential equations and 





boundary conditions to be combined into an integral equation 
for the wave motion; the solution of the integral equation 
provides the intensity in the Far Field. This method could 
practically be used in collecting submarine radiated acoustic 


Signatures. 





ii oLNE hae THEORY 


The analysis of the disturbed region has to be done us- 
ing a different approach than for discrete frequencies and 
one of the most common becomes the correlation technique. 
It is well known what the relationship is between the auto- 
correlation function of random signal and the total average 
power and the D.C. power contained in the signal. Ina 
sense the autocorrelation function 1S a measure of time 
variation, spatial variation and statistical dependence in 
the region of inhomogeneities. 

The study of the statistics of noise gives the possi- 
bility of finding the solution of the problem using the cor- 
relation function instead of a frequency decomposition and 
individual analysis of the frequency components in the ran- 
dom signal. In the process any phase information is lost, 
but the important point becomes the relation between the 
autocorrelation function and the field intensity, which is 
phase independent, and the angular dependence with respect 
to the acoustic axis of the source. 

It is helpful to begin with a non periodical vibration 
problem, with several degrees of freedom, each degree of 
freedom has its own equation of motion and the general equa- 


tion becomes: 





M n = 5 Ge) + BE Gone) - D (vt) COPD 
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where Bae) is the internal force, aCe) 1s the exter- 
nal force, D_(v,t) is the frictional force, and M 1s the 
effective mass of the system. 

Plotting the behavior of one component of x and its 
associated velocity on the Nth phase plane, 1t can be shown 
that the contour is not a closed-loop. Usually it will be a 
complicated non periodic curve confined within some limits 
given by the energy of the system but never repeating itself. 
Of course the N gegree of freedom implies N different frequen- 
cles in such a way that the combination is not periodic. The 
displacement component will not be a periodic function of 
time, so the Fourier Analysis can not be applied. The average 
value of the displacement is also time dependent for short 
periods of time, but if the time is long enough the fluctua- 
tion can be evened out, and the resultant average becomes a 
constant; again it is necessary to use the statistical pro- 
perties of the autocorrelation function, which is defined as 
follows: 


R(t) = Lim [A (Poe se@ey) 8 Ce ret (2.2) 
-T 


T'>00 


The value of R(0) is equal to the average power of the 
Signal over the interval, and this will be the maximum value 
See Cem lie size of “t" Isevery important since the 
rapidity with which R(t) drops to zero as t is increased 
is a measure of the degree of randomness of the oscillation 


of the function. In the case of the spatial autocorrelation 


IE 





function, the correlation distance, (d), follows the same 
behavior. 

The second step will be to study the radiation character- 
istic from a plane surface with randomly vibrating area, the 
goal in this section is to try to obtain a Green's function 
Such that one or more parts of the resulting integral equa- 
tion can be eliminated, and the remaining part easy enough to 
be eeerated. The solution of this problem will be approached 
moe fOour steps: 

Semeciciilation Of the far field due to a point source in 

the presence of a rigid boundary serves to outline the geometry 
and describe the role of the boundary and facilitate the choice 
of the right Green's function for subsequent steps. 

(2) Calculation of the far field due to harmonic motion of 
parts of the surfaces itself which is assumed to have surfaces 
characteristic of a rigid boundary. 

(3) Calculation of the far field due to harmonic motion of 
parts of the surfaces which have a particular acoustic 
impedance. 

(4) Application of results from (3) to randomly vibrating area. 

1. Consider sources in the presence of an infinite bar- 
rier like the arrangement shown in Figure 2.1. If the plane 
is perfectly rigid, the boundary condition is that the par- 
ticle velocity and the glee. Oke the pressure at the surfaces 
is zero, thus, the effect of the boundary plane can be replaced 
by a set of images sources, symmetrically located with respect 


to the plane, both radiating into the region containing Po. 
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The acoustic pressure at @ , defined by the 
vector Fr » for a simple-harmonic source is then: 


1 etkr1 + —2 eikrg),~iwt 


p = -ipck SH ea 


, i (e3) 


where k is the wave number and SH represent the strength 
of the source expressed in the frequency domain. Some 
approximations have to be made in Equation (2.3), assuming 
that Q is at a distance such that the two waves combine 
to form a source cate looks like a single nonsimple source 


= 
located at the origin. The magnitude: of vectors ry and 


= g& 
ry is defined by: 
> => r > 
= sas = _ ee = = ‘e 8 

ry |r aad eC) ae ie haeces 
and 

r, = |r-r't| = r+ 2z cos a-r._ cos ¢ 

2 O O O 


Using these approximations Equation (2.3) becomes: 


ik(r - r_cos ® + z_ cos a) 
~iock S cos (kz. cos a) e 


Poe Dune Ww 


The distance from the source to the plane is an 
important parameter since it determines the constructive or 


destructive interference. For example, if KZ, and kr, are 


14 





much less than 1, the source is much smaller than a 
wavelength and the behavior of the pressure in the far field 
is very similar to that produced by a simple source of 
double strength of the original source. On the other hand, 
alae Zo is equal to a quarter wavelength the separation 
between the source and its image is one half the wavelength 
and its effect is a destructive interference. The same 
effect can be produced if instead of a point source a 
Teele oLe = SOurce exists for which each point has its own 
image. Of course the boundary produces distortion in the 
near field which is not considered in this paper. 

2. Now, considering that the sound waves in the bounded 
medium may be produced by the motion of the plane surface 
itself. Withthe plane rigid at Zz equal to zero, some 
Perron ore the plane may be moving in the 2 direction. 
For a simple harmonic motion ae exp(-iwt) , the 
boundary condition is that the gradient of the pressure is 
ikpcu, (x,y)exp(-iwt) » whereas, for the rest of the plane 
the partial derivative of the pressure with respect to 2z 
is zero. 


The time independent part of the solution of the 


wave equation in terms of the Green's function is given by: 


pyr) = SEBS av, + IRC p(t) - Put o) aa GG] a 
(2.4) 


1? 





Su@eemomere sre ne sources in the medium the volume integral 
is equal to zero and the surface integral involves p and 
the outward pointing normal gradient at z= 0, but the 
feceripe1on of the physical situation has specified only 

the normal gradient. The Green's function should be selected 
in such a way to eliminate the second part of the surface 
integral, keeping only the specified normal gradient of the 


pressure. The Green's function satisfying this requirement 


is given by: 
Tae ier 
cee ee irs ae ee 2 
wr Anr dtr 
O 1 eo 
in this case Pr, = Po > because the source is in the surface 


boundary, thus, the bounded medium Green's function (G) 
is twice as much as the unbounded medium Green's function 


gz the reason is that each moving element of area radiates 


U) 3 
only into the half space z»> 0 , the source and its image 
Mavemcoalesced to form a source of double strength. 
Applying this concept, Equation (2.4) becomes: 


ve) dx dy (2.5) 
-ikoec fFu, (xo »¥5) G 6X5 29 52%, O O 


Py 6 Xo¥ 22) 


= X27 2% 
2ikpc fFu, (x, 2.) by eye 0 dx, dy. 


3. A perfectly rigid boundary is unrealizable. The 


common situation is that the boundary has a particular 


acoustical impedance defined by its density and the velocity 


16 





of sound in the medium. Then, the boundary reacts to an 
incident pressure and every portion of the plane is moving, 
and a different method has to be applied in order to dis- 
tinguish the active regions in the plane, where energy is 
fed into the medium from the passive areas. 

For the passive areas the relationship between p 


and its normal derivative is given by: 


Ce oe | — 
a a wae 5 ey aie Z a C26) 


where 8 is the impedance ratio between the impedance of 
the medium where sound propagates and the impedance of the 
surface boundary. For the Active Areas the normal fluid 
velocity, plus the admittance of the boundary times the 
pressure at z= 0 must differ from zero, this difference 


will be defined as the driving velocity Cu) 


Cu, : oc _ (ikpe OZ a oc ar Bar (2.7) 


To simplify the surface integral, the Green's function 
could be arranged such that the surface integral over the 
passive areas is zero. This is accomplished by a Green's 


f 
function that satisfies: 


dG 
Pertined 


OZ 


= -ik8G G20) 
2 Ww 


aby 





Preweicceraie over the active area is: 
W 
ee dx, dy, Cau! 
Usmmieeede (2.6) am Eq. (2.9) 
9D 
p = - om + ik6p J G ax, qy, G21. 0) 


The expression in parentheses is the driving velocity in 


the active region given by (2.7), thus Eq. (2.10) becomes: 


= —_ Sean 
Py ikoec fu, OF J, ) Swe ‘Ve 5) dx, dy, 2a) 
which has the same form as Equation (2.5), but this time an 
exact representation of the Green's function as was stated 
previously is not possible. However, a good approximation 
valid when the measurement point Q is not closer than a 


half wavelength to the boundary surface and which satisfy 


Equation (2.8) is given by: 


ikr, ikr., 
om cos a! = 3 ]° e @ 12) 
Gy = War, om cos a! Lar. ' 


“where Phemvartables correspond to Fig. 2.1. Although this 
form of the Green's function is derived for the problem of 
ae Oluiemeg@unce ime pne Dresence of @ non=-rigid boundary, it 
may be used since it accomplishes the simplification of the 


problem in hand. [In this approximation there is an image 


18 





at ee » but the wave it radiates has an amplitude which 


depends on the angle a! which the line pr makes with 


2 
the z-axis normal to the boundary. The ray of acoustic 
radiation which reaches the point Q from the source point 

at - by reflection, strikes the reflecting surface at Q 

at an angle at of incidence and reflection. The wave 

which goes directly from Ee to Q without reflection, 

given by the first term, is independent of the inclination 

of ry 3 the reflected wave, given by the second tern, 
depends on the angle of reflection. The factor 8 introduced 
EMemvarlation in the reflected sound due to the boundary. In 
this way the term in per enencscs represents the reflection 


coefficient. ] 


In the far field Eq. (2.12) can be simplified using 


ry = Po: 
ikr, ikr 
_ e COSMO t= ee Cosng' 
aa Tar, [1 + cos a! + g amr, beas at + 5 J (2.13) 


From Eqs. (2.11) and (2.13) the pressure in the far 


field is given by: 


ikr, 


- eeidele CO eh e 
OT Ju (x9 294) cos a! + 8B ry ax, dy, (2.14) 





Py 


Observe that the integral is over the active region 
of the plane which is assumed to be confined to the circle 


of radius "a" about the origin. 


Wy 





In the far field the last equation can be simplified 


neglecting terms of order a/r and smaller and using the 


Seproximation Cos qg' = cos q = Z/r , 
= = cos ct Si i 
ry ig Oe, ¢ + y, sin ¢) sin a 
ae = and tan » x z 
Py te X 


Thus, Equation (2.7) becomes: 


ikr ~ik:W.sin a 


COs.o , e O 
ffu (x sve 


8 + cos a Ane 





dx dy 
(25 5) 


2 ~  -ikpe 2 


where the vector Wo lies in the xy plane, with Xy > ve 


as component; and vector k has the r direction whose 


eemoaonences are: 


Ky = k cos 9 Sina 
Ky = aig oi) Shihal ve 
kK, = k cos a 


A close study of the integral part of Equation (C2 Sy) 
shows that the angular distribution of the Far Field is 
proportional to the two dimensional Fourier Transform of the 
distribution of the driving velocity in the active region, 


which is defined by: 


ZAD) 





xO 
a Ie ip Ue ay) € _- lca 


bare O 


1 00 -k x “KY 6 


U (kak) O 


—O 


LK x otk, 
Ie U lk, sky) = dk, dk. 


= OO 


Be 2) 6) 


o 
£ (a5) . 8m cos a 


Bt cos a Upik cos ¢ sina, k sin @ sina) (2.16) 


Thus Equation (2.15) takes the form: 


VG Te 


Cpe ef a4) | Cam 


The angular dependence is given by the factor i A eelal(elal 
measures the relative amplitude of the pressure with respect 
to the center of the active region. 

Several properties of the Far Field can be studied. 


If the coordinate axes are rotated to Xt such Chay Xo 


0? 
axis is the intersection of the rz plane, the appropriate 


transformation is given by: 


m4 
ii 


a X cos 6 - YG Sin 9 


Nie Xo sin ¢ + to Se d 


Under this transformation the time independent particle 


velocity becomes: 


oe 





1 “kX, Sin a 
U Ck cos ¢sin a,k sin ¢:sin a) = Paice é dx. f U(X) ou 


This equation shows the one dimensional Fourier 


transformation in the X, direction, of the average: 
u,(w,X_) = — ee oa dy 
ye fa! OT ft) "= ie O 


of the particle velocity U6 Xo 2%) in the x, direction. 

It can also be observed that Uy (wsX,) iS waehumeb on 
‘opm X, » the distance in the direction of the measured 
meant @ ; it is also a funetion of 6 , the direction of 
the rz plane, and of the frequency w/2m7 . 

If the particle velocity u,(x,y,t) is an arbitrary 
function of space and time and u (x,y) is its time dimension 
Fourier transform, then the pressure field at a point Q (r,a,9) 


is the inverse transform of ar given by Equation (2.17) 


S -iw(t - =) 
p(r,t) = = f f Caso) e du) (ris) 


Now, applying the rotation of axes as defined before, 


Equation (2.18) can be expressed as: 


= BVM ee) © i(-)X_ sina 
cs wt. COSO% =f Ly ° dw f aa ax 
p(r,t) = omr 8 + cosa Coane . = O 
ef We -- rp eM? wx ,Y ,Z) a 
_ ye: ts (2.19) 
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- £ , cosa 9 + _i 
Get) = [lf aX Say, u(X.¥.t + ox sina =) 


8+ecosa ot 
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This equation shows clearly some of the properties of the 
Far Field. An appropriate summation of the normal acceleration 
of the active area gives the total field at the point of 
emcverest. Also the contribution to the pressure at r, at 
meme t , from the line Se > “m7 < — < © , depends on 
the normal acceleration on the line at time 
l/c) X sin(a)-r/fe , at time r/c - (X/c) sin(a) earlier, 
which is the time required for the sound to travel from the 
mene to point Q. 

4h, All the formulae developed so far can be used to 
compute the mean intensity at r,a,?, produced by an area A 
of the boundary which is vibrating with a random distribution 
of velocity. The particle velocity in the active area A is 
a function of position as well as of time. In this regard, 
it becomes a random variable and a way to describe its 
BeWavirormrs Using ifs statistteal properties. In this 
particular case the autocorrelation function is used as a 
means to describe the behavior of the random source. 
Considering the particle velocity as a function of position 
and time in the active area, the autocorrelation function 


is defined by: 
RCEsn,t) = = ffS us(x,y,t) u,(xté,ytn, ttt) dx dy dt 
where the integration is over the active area A , and T 


is over a long interval of time. Also the Fourier transform 


of the velocity is given by: 


eS 





1 twt~ik, x-ik 
U,, (Kok) = ie ff[fu,(x,y,t) e dx dy dt 


From this equation the space time spectrum level of 


the particle velocity is: 


” tuwt-ik_&-ik_n 
i os + IPP Js S Syncs WE Z Y de dandt (2220) 


| U 
64n 


itiecae far field 16 is Known that the intensity is 
eyen by the real part of the pressure times the particle 
velocity, but also the velocity is related to pressure for 


Piemrelation u = p/oc , thus using Eq. (2.16) and Eq. (2.17): 


Dy ae. 2 : 
ee) Ee A he Anvk' cos” a |U (k sinacos ¢, k sinasin ) | 
pc ne Ie ae Gee a} - 


C2521.) 


This equation is related to the autocorrelation 
function by Equation (2.20). The problem is now reduced to 
getting the form of R(&,n,1t) in space and time. A 
reasonable assumption could be a Gaussian waveform where the 
value of R falls off smoothly as its arguments are increased. 


Assume that the autocorrelation is given by: 
~—5 Cae iy - = (w,t)* 
ie ed 


BB olig ee) |U e C252) 


O 
where U.S is the mean square particle velocity of the active 
area averaged over time interval T ; d stands for the 


correlation distance along the plane, WO, is the upper 


a4 





angular frequency limit of the spectrum, which also can 
give the correlation time, and the lower limit of the 
wavelength. 

Using this value of the autocorrelation function 
and equations (2.20) and (2.21), the intensity field was 
evaluated: 
gg sre lion? a REF auto KE (2.29 


I(w) = pe -=> > 
7 ré ew. |B + cos a| 


BUG vdeoes onffy fer a specifie frequency; in order 
to get the total field intensity it would be necessary to 
eonsider the contribution of ail the frequency components 
of the random signal. <A logical way to accomplish this is 
to integrate Equation (2.23) with respect to frequency. 
Note that the integrand has a maximum when the exponent 


is zero at: 





me=ot’] jG & “a> sin’ «)*/* 
ig n 
a 4. Wal cosa = (ik, 
ae . Oo 
Lr) = ar ff T(o)dw = pe 5° ——————3 D2 aa2 
a emr =Ss«d{B + cos a| [1 + (k,d)“sin® a] 


(2.24) 


A simple analysis of this formula shows that the 


intensity is proportional to the impedance of the mediun, 
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Piicomememmecaesduare particle velocity in the active area, 
times the ratio between the active area, A , and the area 
Seoede Ralf sphere of radius r, times an angle distribution 
function , which depends on the maximum frequency of the 
spectrum (inverse of the correlation time), the speed of 
sound and the correlation distance, d. 

Trying to apply this theory to practical cases, two 
main problems arose. The first is that in real situations 
the frequency spectrum is band limited, usually due to the 
physical characteristics of the source. Under this 
consideration the limit of integration should have an upper 
bound (UL) and a lower limit (LL), in this form, Equation 


(2.23) becomes: 


Iu. |2cos2 UL -3)*sinXa- 5 (4)° du 
AT 4-5/2 @! Se ae 2 : 
b G 2 2 m pe A 

r cw. |B cos a LL 


C22 >) 


There is no analytical solution for this equation, so 
a numerical method of integration had to be used in order to 
get the intensity in the Far Field. 

The second problem is the difficulty in measuring the 
autocorrelation function of the randomly vibrating area, 
which is related to the particle velocity in the active area. 
Instead it can be assumed that if the area close to the 
active region is vibrating with similar characteristics, 


then the impedance ratio becomes unity and the autocorrelation 
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Mmaoemicanwevte la 2ivyes the mean square particle velocity; 


Equation (2.25) becomes: 


Din ee um - aC) sin® a - HB)? au 
A 3/2 4 JU,|~ cos” a ; n 
2) = 2c 5 (21) ae we 
r cw {1 + cos 0. | rL 
(2.25) 


This equation is angle dependent, thus the intensity 
in the Far Field depends on the angular position with respect 
to the acoustical axis. In this way it is possible to 
evaluate the beam pattern. Observe also that WO, should 
lay inside of the limit of integration or at least equal to 
the upper limit, otherwise the process would show a 
discontinuity which is not possible because the frequencies 
eemermimuoine to the field are given by the limit of 


integration. 
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III. EXPERIMENTAL RESULTS USING THE PLANE MODEL 


To validate the theoretical model described in Section 
II, Far Field and Near Field measurements were made on a 
horn speaker in an anechoic chamber. The actual beam pattern 
in the far field of the horn-speaker was taken and also the 
autocorrelation function in different points of the near 
field were obtained in order to get the mean square particle 
velocity in the Near-Field and be able to predict the Far 


Field Beam Pattern. 


A. FAR-FIELD MEASUREMENT 

The normalized Far Field beam pattern was obtained using 
the experimental set up shown in Figure 3.1. Since the 
source was driven by random noise the distance of the near 
field is variable. The criterion for near field is usually 
given for one particular frequency, as the frequency is made 
lower and lower the near field also decreases and eventually 
disappears for practical purposes. That is the reason for 
the presence of the high pass filter which limits the low 
frequency component to 1300 Hz in order to get an appreciable 
near field region. This 1s necessary for obtaining far field 
and near field measurements under the same frequency spectrum. 
The signal is then amplified and applied to the horn speaker. 
The microphone was placed in the far field at a distance of 


9-8 meters from the source, the received signal was amplified 
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and the output connected to a true RMS Voltmeter which also 
generates a DC voltage output proportional to the deflection 
of the needle for the selected scale, this voltage was ap- 
plied to the Y axis of the plotter. The xX axis of the 
plotter received information from the position indicator as 
a DC voltage proportional to the angular displacement of the 
acoustic axis of the horn speaker. The normalized pressure 
amplitude versus angular position is shown in Figure 3.2. 
This was the actual Far Field Beam Pattern which had to 
be compared with the beam pattern evaluated using the mathe- 


matical model described in the previous section. 


Bee NEAR-FIELD MEASUREMENT 

The intensity in the Far Field, as defined by Equation 
(2.26) depends, among other factors, upon the mean square 
particle velocity. The mean square particle velocity is a 
source quantity which was shown to be related to the auto- 
correlation function. The evaluation of the factor is 
possible by using an array of microphones over the entire ac- 
tive region and the summation of the received signal fed 
into the correlator. The microphone array was not available, 
instead a single microphone was used to take a measurement 
each 5 cm in the vertical and horizontal direction trying to 
cover the active area as shown in Figure 3.3. The autocorre- 
lation values obtained at each single point is then averaged 
and the resultant value was used to describe the mean square 


particle velocity. 


30 





*daxeeds udgoy *ugezzed weegd PTET UeJ poeunseoay Z°€ sansty 


1 


/ t 





bs ey ee ee ee ee ee ee - 2. - - . = - > = ~ 22 q-- « - ~ 
5 - 


[sboibopjx 
A = eee eas ae 


: ; { ‘ 
) a 


; | 5 ’ 
\ t 
“ a 


. i 5 : 
2 re = - as en i oil ) _— a et erent 


1 . 
Re rene pees 


= 

Bs esl a eee 
' ; 
: { 


- - —d4 -_—_ +0 
) 


oe ye oe 


! 


’ ' 
as we b ee ee ene °°. 


a . ‘ . 
—— “eee oe en ee SE Cate ee ——— pe - A eg enn * 
’ : 
: ' 


| 
| 
4 
| 
| 
if 
| 
je om 





; ] 
es oe tie af eee fe = Be Semen 
. 
: / . ' apna TTduy pez tTeuon : ae 





OS 


ou 





LJ 
a 
= 
Lt 
© 
2 8 
= 
= 
© 
= 


HORN SPEAKER 


SYA 


BOUNDARY 





Near field schematic representation of 


the position of the microphone. (View facing horn) 


Figure 3.3 





The correlator was a SAICOR model SAI-43A, which was 
used in the autocorrelation mode, the function is determined 
Simultaneously at 400 incremental lag points, the completed 
correlation function is displayed at one time, the function 
is provided with incremental lag or time delay value rang- 
ing from 0.2 microseconds to 1 sec resulting in total time 
delays from 80 microseconds to 400 sec. Precomputation de- 
lay of 800 lag values selected in 200 lag increments allows 
the autocorrelation function to be viewed symmetrically 
about zero or up to 800 lag values removed from zero. The 
averaging is accomplished digitally in a linear fashion with 


fixed summations ranging from ae EO pel 


in binary steps. 
The quantization is 8 bits per channel. 

The SAI-43 operates as follows: one input signal is 
applied to either channel, and internally routed to the 
second channel, the signal is sampled at the rate indicated 
by the operator-selected sample increment and then is quan- 
tized in both channels. Channel I is the channel which con- 
tains the 400 word memory. Channel II contains the one word 
memory - the present sample. At the time called "Present 
Time", the correlator performs the multiplication of the 
present value word of Channel II by each of the 400 words 
in Channel I. The product "0" times "0" becomes an entry in 
R(0), the first point of the autocorrelation function the 
product "0" times "1" becomes an entry for R(1) and finally 
"0" times "399" becomes an entry for R(399). At the end of 
this arithmetic cycle the correlator has performed one calcu- 


Terentia aeheso-called bin of the autocorrelation function. 


a1. 








The product in R(0Q) is the first entry for the zero-lag auto- 
Gorrelation value, while the product in R(399) is the first 
entry for the correlation value at a lag in Channel I equal 
to 399t. These first entries, resulting from one sample in 
Channel II are expanded to the full correlation function by 
taking repeatedly updated value of the sample in Channel II 
and correspondingly updated values in Channel I and then 
averaged. In effect the 400 word memory slides ahead one 
word and the arithmetic cycle begins again. For all t incre- 
ments that are equal or longer than 500 usec the SAI-43A is 
Said to operate in real time. But for increments less than 
500 usec the figure is different. Assume an increment less 
than 50 msec, it could have completed eight cycles in 400 
msec instead of one, therefore it 1s eight times slower than 
the real time instrument. 

The output of the correlator can be displayed in an 
oscilloscope or in an x~-y plotter, the voltage output in the 
autocorrelation mode can be determined by the following 


formula: 


Sg) = x (VOLT) dBA/10 


5 eet 


Coy 


where x 1s the voltage as ee as the output of the SAI-43A 
dBA is the dB setting of the input channel. R _ is the 
correlation value in (volt)? to be determined. 

For illustration purposes Figure 3.4 and Figure 3.5 show 
Piowaticocouretlatlon funetion, R(0,0,0), for the point of 


maximum intensity, (center of the speaker), and the point of 
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minimum intensity, which correspond to the corner of the 
boundary of the active region, the relative power between 
them is more than 20 dB. | 

The physical arrangement in the anechoic chamber and the 
instrumentation for taking the autocorrelation function at 
each point in the near field is shown in Figure 3.6. The 
autocorrelation of the band limited noise becomes the sin x/ 
x, function, the amplitude of the function is variable depend- 
ing on the position of the microphone. As described before 
it is also possible to get the power at each point of the 
autocorrelation function. 

The criteria used for determining the active area was on 
the basis of the differences in power between the center of 
the horn up to the point which had 20 dB less power, which 
corresponds to the boundary of the active area shown in 
peepure 3.3. 

The other important measurement was the frequency spectrum 
of the received Sono le wires ven an FLeure 3.7. This 
was obtained using the basic setup but with a wave analyzer 
instead of the correlator. The frequency sweep was made very 
Slow in order to get a good average. The main use of this 
Figure was to provide the proper limits of integration for 
Equation (2.26). The limit of integration became a quite 
important parameter since it modifies greatly the form of the 
normalized beam pattern; thinking in terms of a transducer, 
it is known that higher frequencies make the beam narrower 


and the low frequencies have a wider beam, so looking at 
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Figure 3.2, the actual beam pattern, it is logical to assume 
ee the main contributions are given by the low frequencies, 
and that is exactly what the Frequency Spectrum shows in 
Figure 3.7. The lower limit of integration is assumed to be 
1100 Hz and the upper limit is 5000 Hz, the contributions of 
the other frequencies was small; this also showed that the 
Speaker has a low frequency response, although it covers the 
important band of the human hearing. 

Another critical factor in the Near-Field measurements 
was the correlation distance, since its effect in the width 
of the beam is great. The symmetry of the source caused some 
difficulties; it was found that all points located at the 
same distance from the axis of the horn had a very good cor- 
relation, regardless of the distance between them, but along 
the same angular position a distance of 12 cm was enough to 
diminish greatly the spatial autocorrelation function. This 
Suggested that the correlation distance for this experiment 
is a vectorial parameter, but was used like a scalar along 
this work. It 1s useful to point out that in a real situa- 
tion with an asymmetric source driven by several sources with 
different frequency bands, like a submarine in which each 
section is affected for different set of frequencies and also 
the physical configuration changing in space, the autocorre- 
lation distance has to be easy to measure experimentally and 
its vectorial characteristic should disappear leading to a 
Scalar pa nenerer: Otherwise the pattern must be built up 


from several sources. 
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Along this experiment a correlation distance of 12 am 


was used. 


C. RESULTS FROM NEAR-FIELD MEASUREMENT 

The first step in the data processing was to get the 
average mean square particle velocity. Using the data col- 
lected in the experiment and the formula given by Equation 
(3.1) it was possible to evaluate the power at each point of 
the active area, this power was averaged, the entire process 
was done using an IBM-360 computer program in fortran IV 
language. The value obtained in this program represents the 
mean square of the particle velocity expressed by the 
‘Equation (2.22), where the value of the autocorrelation 
funetion 1s at zero spatial and time shift. 

The approximation used for that evaluation was, of course, 
an additional source of inaccuracy in the experiment, but 
should not affect toomuch the final result, because what is 
really affected is the intensity field. The normalized beam 
pattern is practically unaffected since the relative power 
between each angular position was the subject of interest. 
In the case of trying to obtain a radiation pattern, which 
1s beyond the scope of this work, then Uo? becomes a very 
important factor since it affects the Power Level of the 
Signal. 

The solution of Equation (2.26) is the actual prediction 
of the Far Field Beam Pattern. Again the computer system 
was used for performing numerical integration, a library 


Sub-routine called DQATR was used. It was necessary to use 


41 





double precision in order to obtain a reasonable accuracy. 
The auincaopessige uses the trapezoidal rule in connection with 
Rombergs principles for evaluating the approximation of the 
iieesral value of the function. 

The resulting computed Far-Field is shown in Figure 3.8 
in comparison with the actual Far-Field. Note that the 
agreement between the predicted and the measured beam pattern 
is fair in shape but slightly out in values, the actual beam 
pattern shows a wider beam that the predicted one. Several 
factors can produce that difference, the correlation distance 
could be one of the main factors due to its specially vari- 
able characteristic in this experiment, as it was explained 
before. Also the limits of integration assumed an important 
role, strong high frequency components could produce that 
narrow beam. 

In order to study the variation of the function with 
respect to the limits of integration more calculations were 
made reducing the upper limit to 4250 Hz, all other parameters 
were kept constant. The result of this calculation is shown 
in Figure 3.9. This time a better agreement was obtained 
between the actual and the predicted beam patterns. A wider 
beam was obtained but still showed a small difference with 
the actual beam pattern. 

The effect of the correlation distance is shown in Figure 
3.10. The limits of integration were kept at 4250 Hz and 
1100 Hz, and the correlation distance was increased to 15 cm 


Pmeiemeahemlattonmen it Gan be seen that @ small variation in 
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this parameter affects greatly the predicted width of the 
beam. This Figure also shows that the beam pattern is more 
sensitive to the correlation distance than to the limits of 


‘integration. 
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IV. CONCLUSION 


The possibility of predicting the far field beam pattern 
from randomly vibrating surfaces has been demonstrated. It 
1s feasible to obtain the intensity field and the radiation 
pattern from near field measurements of distributed sources. 
Although the theoretical model developed in Section II was 
tested experimentally with satisfactory results, more work 
needs to be done in both areas. The amount by which the beam 
pattern depends on each parameter makes it difficult to state 
a solid argument for describing the range of validity of each 
one. 

Hieweconpuredmbeam paiteriewas Obtained Solving equation 
(2.26) using a computer program for solving the integral by 
numerical methods. The limits of integration were chosen by 
means of the Frequency Spectrum, which was taken using a wave 
analyzer with automatic frequency sweep. The maximum range 
considered was from 1100 Hz to 5000 Hz. It was assumed that 
the frequencies producing the field were inside of this fre- 
quency band. 

The mean square of the particle velocity, which has 
energy information, was shown to be related to the autocorre- 
lation function. A correlator was used for that measurement 
and the power contained in the signal was easily calculated 


with equation (3.1). Note that in this case a correlator was 
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Nscasputeaetrue RMS yoltmeter can also be used to measure 
the power of the signal. 

The correlation distance was also measured with the 
correlator, a pair of microphones were located in the active 
area and the two dimensional (spatial and time) correlation 
function was observed. A correlation distance of 12 cm was 
used for calculation, but the measured value was found to 
vary with location of the microphones. The calculated beam 
pattern is very sensitive to this parameter. Its effects 
are bigger than those due to the limits of integration. 

Different configuration of sources should be tested ex- 
perimentally. Sources should be excited for different fre- 
quency bands in such a way as to produce an active region 
With a well defined correlation distance. This will solve 
one of the main problems of this experiment. Also it is 
Suggested that the paegemnes spectrum of the configuration 
should not necessarily be continuous, as in this form the 
arrangement will be closer to real cases. 

The theoretical model does not consider the spectrum 
level of the signal for the different frequencies. This 
should be tested more carefully in order to state a rule for 
the limits of integration. In the case of discontinuities 
in the Frequency Spectrum the model fails to give a straight- 
forward solution, although it looks reasonable to perform a 
Summation of the contribution of each frequency band. cee 
was not possible to verify this concept with the experimental 


setup. 
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Agreement was generally good between measured and com- 
puted patterns, especially when the upper limit of integra- 
milenewas reduced to 4250 Hz. The actual beam pattern was 
slightly broader than the computed pattern in all cases. 
This is not surprising since the far field beam pattern de- 
pends on the frequency response of the source and the limits 
of integration were chosen in a rather empirical way, so the 
computed beam pattern could include a bigger frequency band 
than the real one and the width of the beam would be too 
narrow. 

The technique developed here might also be useful in 
measuring the radiated noise pattern of submarines. Some 
additional considerations should be taken into account, and 
the theoretical model should be validated by actual measure- 


ment in anechoic tanks. 
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APPENDIX A 


COMPUTER PROGRAM FOR EVALUATING THE FAR FIELD BEAM PATTERN 
USING THE THEORETICAL MODEL DEVELOPED IN SECTION II 


PROGRAM TO EVALUATE THE BEAM PATTERN OF A HORN 
SPEAKER DRIVED WITH RANDOM NOISE. 

FOR PROGRAMING PURPOSES THE FOLLOWING SYMBOLS HAVE BEEN 
[Of DPS 

Ci — CHARACTERS TEC MPEDANCE OF THE MEDIUM 

be= CORRELA LICN pis TANCE. 

Bata = ANGLE Wilh RESPECT TO THE ACOUSTICAL AXIS. 
T = AVERAGE OF THE MEAN SQUARE OF THE PARTICLE 
VEC ar tb SOLVE hEGlLON: 

XL = LOWER LIMIT OF INTEGRATION. 

WN=UPPER LIMIT OF INTEGRATION. 

A AREAS OPS the ACTIVE REGION. 

C Seo ee or SOUND IN THE ALR. 


CALCULATE THE INTENSITY IN THE FAR FLELD 


C2'Ca C2 Go") (9 9) 9) a 9 ON ee eae 


IMPLICIT REAL#8 (F- H ,Q-Z) 
REAL LABEL/' 7, 

REAL*8 TITLE (12) 

REA. cme ee Coeo) Dc .cChoCS 
DIMENSION P(73),BETA(73) 
EXTERNAL F 

COMMON D,C,RBETA,SN,PI 

PI = 3.1415926 


D = O0.t2 

EPS = 0.001 

WN = 2.*PI*5000. 
NDIM = 999 

A =0.5**2 

T = 0.062 

CI = 415. 

C = 343. 

R= 5.8 

NEE SPIER) (OV 
Si> CLaA* 1 *(DA*2 ) 
WP = 4250. *2. * PI 
S2 = (R**2) *(C*#*2) *wWP*((2.*PI) ¥**#1.5) 


DOO Tis. 1 237 

BETA = BETA(I)*PI/ 180. 

SN = (DSIN(RBETA) ) ##2 

CS = DCOS(RBETA) 

CALL DQATR(XL,WN,EPS,NDIM,F,Y,IER, AUX) 
P(I)=(SI*(0S#*2)#Y) / (S2*((1.+08) **2)) 
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Cera) OG 


BETA(74-I) = -BETA(I) 
P(74-I) = P(I) 
WRITE(6,100) Y,IER,BETA(I),P(I) 
10 CONTINUE 
100 FORMAT(/'Y =',E16.7,5X,'IER=',15,5X,'BETA=',F7.2,'INT= 
1 SeEt6 7) 


CALCULATE THE NORMALIZED AMPLITUDE OF THE PRESSURE FAR 
FIELD BEAM PATTERN. 


Z = P(37) 
DO 20 I=1,73 
ROS «CRG 72) 2075) 
WRITE(6,200) BETA(I) ,P(I) 
P(I) = P(I)* 2. 

20 CONTINUE 

200 FORMAT(/,'BETA =!,F7.2,13X,'NORMALIZED AMPLITUDE =',F 
WS) 
READ(5,50) TITLE 

50 FORMAT(6A8) 
Chili PRoviGieeRe Ee ENO On LAB TlmiEn.~20.,0.5,0,4 0,2, 
WEN Is 0) 5 IySGE ) 
WRITE(6,400) LAST 

400 FORMAT(/'tLAST =',I2) 
STOP 
END 


FUNCTION F(W) 
IMPLICIT REAL*8 (A-P,0-Z) 
COMMON D,C,RBETA,SN,PI 


BLIP = O75 
BFP = 4250. 
ERP = 24 


WP = BFP * PI * BRP 

F =(W*#2)/((DEXP(BIP*((W*D/C)**2)*SN))*(DEXP(BIP* 
1((W/WP) **2)))) 

RETURN 

END 


as 
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